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ABSTRACT 


An increasing emphasis has reeently been plaeed with fighting non-national and 
irregular forees. Single entity attaeks on shipping and transportation units have replaeed 
attaeks on fortified or established military positions. The supply ehain from sea to shore 
has beeome the target of opportunity to disrupt the Global War on Terrorism and put 
further lives in danger. Autonomous eontainers will play an essential role to deliver 
logistieal supplies to waterborne littoral vessels in order to maintain station and eomplete 
military operations all while eliminating the threat to human life as the eontainers will be 
programmed to deliver supplies to a specified loeal in a reasonable timetable; vessels 
sueh as Riverine Warfare patrol eraft, Speeial Operations eraft and Coast Guard seareh 
and rescue boats and their erews. The researeh to be eondueted will foeus on the 
seakeeping eharaeteristics of an autonomous sustainment eargo eontainer and the 
feasibility of its deployment. Established geometrie data will be used along with 
ehanging loading eharacteristies and ride effects. The in depth analysis will be foeused 
on the responses of the eontainer in varying sea eonditions and at varying loads to see if 
further refinement of the design or polieies concerning loading and deployment may be 
required. 
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1. THE AUTONOMOUS SUSTAINMENT CONTAINER 

CONCEPT 

A. BACKGROUND 

Throughout the history of sea to shore engagements, supplying goods from units 
at sea to forces ashore has met with several shortcomings. Most notably, in the World 
Wars, Korea, and the recent conflicts in the Middle East, continuing to support troops 
ashore from forces afloat carries with it high risk, both to materiel and personnel. The 
current combat environment has placed further strain on this link of the supply chain as 
non-national forces look for easier and faster ways of harming the war effort than 
attacking fortified and/or heavily manned emplacements and convoys. A new method for 
delivering the necessary sustainment goods to forces ashore without endangering lives or 
providing easy targets of opportunity for these opposing forces needs be found. 

B. HISTORY OF SEA TO SHORE SUSTAINMENT 

I. Traditional Methods 

To provide supplies to forces on shore, the method used throughout history 
required the supply ship to dock with a physical structure designed to move goods from 
ship to shore. As ships, especially cargo ships, grew in size, so did the required dock 
size. The dawn of the 20**' century saw cargo ships that were too large to traditionally 
dock that moved supplies from anchor via smaller utility boats. This method still 
required a ship to berth at a physical dock to unload supplies. World War II saw the 
advent of supply sustainment without the physical dock. Army ships simply pulled up to 
the beach and offloaded supplies. This method of beach landing to resupply troops on the 
beach has remained the preferred mode since that time. As seen in Figure 1, U.S. troops 
deploy from an Army LCV on D-Day 1944 compared to U.S. and Royal troop landing in 
2007. Although the commissioning dates of these two vessels are significantly different, 
the methodology remains primarily the same. The same is true about the delivery of parts 
and supplies from sea to shore. 
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Figure 1. U.S. Marine landing at Omaha Beach, 1944 (left) (From: 10) Joint 

U.S. and Royal Army troop landing 2007 (right) (From: 11) 


2. Non-traditional Techniques 

There are points in history, where eireumstanees prevented the use of tradition 
supply delivery methods. Winds, waves, prevailing eurrents and even single individuals 
were used to eireumvent the need for both berth and beaeh, or even for the ship to eome 
in eontaet with land at all. “...the Tokyo Express was [having diffieulty delivering 
supplies]...so the measures now being taken were desperate ones. Destroyers tried to 
supply these Imperial troops by making high-speed runs to Guadaleanal and dropping off 
strings of lashed-together drums into whieh supplies had been sealed. Barges from the 
island then were to tow these drums ashore. The troops managed to retrieve only about 30 
per eent of these supplies that Tanaka's destroyers east upon the water. Tanaka's 
destroyers raeed down The Slot on 3 Deeember and dropped strings of 1,500 drums. But 
the island troops managed to haul in only about 300 of these from the waters off 
Tassafaronga and Segilau. On 11 Deeember, Tanaka himself led one in his speedy 
Teruzuki... and dropped some 1,200 drums of supplies off Cape Esperanee and then 
headed north again. Only some of the drums were reeovered by the troops ashore...”[5] 
The plight of the Japanese troops eontinued until the islands were finally conquered by 
U.S. troops. Eurther stories of the Tokyo Express include using divers to clandestinely 
recover supplies floated by Japanese destroyers and midnight barge runs through U.S. 
patrolled waters. When the traditional ways begin to fail, it is essential to make every 
effort to solve the problem with something new. 
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3. Unmanned Ships in History 

There have been a few instances of unmanned and semi-autonomous vessels 
taking part in important naval operations where manned vessels simply would not have 
been as capable. One of the most famous uses of unmanned craft occurred during the 
battle between the Spanish Armada and the Queen’s Fleet off of Calais. Sir Francis 
Drake and John Hawkins hoped to pull the Armada out of safe harbor using the then 
famous hellburners of Antwerp, ships that had been made into tremendous bombs, as the 
tool in which to bring the Spanish fleet into combat. “Drake offered a ship of his own, 
the Thomas of Plymouth, of two hundred tons, and Hawkins one of his; and with the 
mounting enthusiasm six more were recruited...the captains scattered to the work of 
getting them ready and stuffing them with everything available to make a blaze.”[3, 
p271] “But not only were spars and sails and rigging left, since the ships would bear 
down on the anchorage...under full sail, but all the ships’ guns were left double-shotted 
to go off when the fire made them hot enough, to add to the terror if not the destruction of 
the enemy.”[3, p271] “.. .the Spanish anchorage could see plainly eight tall ships with all 
sails set and lines of fire beginning to run up their rigging, driving straight towards them 
with wind and tide.”[3, p277] “...the line was driving along so fast, pushed by a strong 
wind, a spring tide and the set of the Channel current all working together...”[3, p277] 
“...the double shotted guns heated now white hot, began to go off, spraying their shot at 
random over the water and by the force of their recoil sending up a fountain of sparks to 
be blown among the boats...the six remaining ships swept past and bore down on the 
anchored fleet, the sound of their exploding guns heard above the roar of the flames and 
the fountains of sparks shooting skyward.”[3, p277] “ ...most captains simply cut their 
cables and ran before the wind, scattering, some here, some there, as if they were as much 
afraid of one another as they were of the fire-ships. The strong set of the current and the 
rising gale swept the whole disorderly mob out through the straits and on toward the 
sands of the Flemish coast.”[3, p278] The autonomous fire-ships did little physical 
damage to the Spanish fleet, they were after all not really the hellburners the Spanish had 
feared, but when dawn approached the Armada was in complete disarray. Unable to 
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reform, they were harried for days by the Queen’s fleet until able to flee northward. This 
marked the end of the Spanish Order and set about the beginning of the English 
domination of the sea. 


4. The Emphasis on Unmanned Craft 

With the development of new eoneepts in military operations and reduetions in 
manpower of our military forees, the promotion of autonomous systems has been pushed 
to the forefront. Reeent improvements in wireless and satellite eommunieations, sensor 
and data relay systems and remote operability has made unmanned vehicles an available 
option for easing the manpower burden felt by much of the military. An autonomous 
vehicle is a self-piloted vehicle that does not require an operator to navigate and 
accomplish its tasks. As the necessity to protect human life grows, the need for 
autonomous vehicles will grow as well. The U.S. Air force has proven, in combat, the 
range of uses for unmanned vehicles. The Air force Predator has run countless sorties in 
Afghanistan and Iraq with missions ranging from simple surveillance to offensive 
surgical strikes. The pilot does all this from a remote operating station well away from 
the threat region. Building off of the advances and successes of the Predator, Global 
Hawk will fulfill mission roles from long range surveillance and tracking to bombing 
strikes against hardened, buried and bunkered targets. Further UAV systems are in the 
works including the Boeing stealth UAV, the MQ-8B Fire Scout helicopter UAV, and the 
autonomous programmable Silver Fox and Scan Eagle UAVs developed at the Naval 
Postgraduate School. These advances have spread to surface and subsurface craft as 
well. Current systems in development include the BAE system Protector USV for Coast 
Guard applications, and the autonomous Sea Fox USV and Remus UUV. 


C. THE DIFFICULTIES OF THE SEA-SHORE INTERFACE 

The sea shore interface is an area rife with hazards. Shallow and shoal waters are 
home to shipping dangers such as reefs, mine hazard areas, capsizing surf conditions, and 
underwater rock outcroppings. The shoals are further threatened during combat as the 
threat level increases with enemy activities. For the Marines at Guadalcanal the problem 

was service support facilities. “Without port facilities, supplies reached the troops only 
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after a series of time-eonsuming and labor-intensive equipment transfers. Supplies were 
first unloaded from Navy ships offshore into lighters for the trip to the beaeh. There 
Ameriean serviee support personnel transferred the tonnage to trueks that hauled it inland 
to several dumps on roads under construetion. From the dumps supplies had to be hand 
earned, by both Amerieans and native laborers... As the fighting moved further inland 
the distanee between dumps and front line lengthened and road building could not 
progress as fast as assault units advanced... American troops temporarily solved the 
distribution problem by using the many streams and rivers on the island. Loading 
supplies into small boats...Americans pushed the craft through the water as close to the 
tactical units as possible.”[4, p26-27] The L* Marine Division faced problems from both 
the shoals and the enemy as they tried to get supplies onto the beach at Inchon. “At 
Inchon, the lack of adequate port facilities made it essential that the naval logistics forces 
rapidly move supplies ‘over the beach.”’[6, p40] While the UDT cleared obstacles and 
the Navy Seabees built a large pontoon dock and reconstructed port facilities, eight LSTs 
beached themselves, under fire, and began the task of disgorging tons of supplies and 
munitions. “After the evening landing, some of [Captain Watson] Singer’s sailors and 
[General Henry] Crowe’s marines had to unload, organize and distribute the supplies 
from the eight LSTs temporarily immobilized on Red Beach...”[6, p40] The precarious 
condition of these beached vessels can be seen in Figure 2, but the surf and sand wasn’t 
the only threat to these vessels or their crews. “Enemy snipers in Inchon shot at Crowe’s 
leathernecks and Singer’s sailors, outlined as they were in the glare of the floodlights 
installed by the Seabees. Nonetheless, the men continued to toil throughout the long 
night to accomplish their vital mission. All of the LST captains withdrew their emptied 
ships on the morning’s high tide.”[6, p40] 
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Ph«U) # 80-(i-42(K)27 I.STs unloading ai Inchon. 15 Scplcmbor 1950 



Figure 2. Four “high-and -dry” LSTs disgorge supplies onto Red Beach on or 

about 16 September (1950) (From: 6) 

D. THE ASCC MISSION 
I. Purpose 

Waterborne logisties delivery systems sueh as the meehanized Landing Craft 

Utility (LCU), Landing Craft Air Cushion (LCAC), and High Speed Vessels (HSV) have 

been used for logistieal delivery for many deeades. All three of those types of delivery 

vehieles are very eapable at providing logistic support with speed and efficiency; 

however, they all provide the opportunity for the loss of human life. The ASCC is 

designed to complete the same missions with equal or better speed, cost and efficiency, to 

effectively remove the need to task those particular vessels with logistic support 

missions. The overall main purpose of the ASCC is to replace the need for the logistical 

interface of sea basing. With the reduction of sea basing comes with it the reduction of 

bulk cargo constraints and cargo breakouts, manpower and the addition in time savings. 

The mission application has immense appeal to Joint Logistics over the Shore (JLOTS). 

The ASCC goes further; to reduce the threat to crew and service support personnel 

aboard both the sea based logistics ship and the delivery vessel itself The reduction in 

manpower requirements, enhancements of personnel safety, and reduction of local 

infrastructure requirements are its greatest attractors. 
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2. Mission Profile 

Figure 3 shows a hypothetical coastline where troops positioned ashore are 
expecting logistic support. The opposing forces have mined the coastline and set up 
surface to surface missile sights overlooking the bay. Under conventional sea-shore 
support missions, the logistics craft would have to close to within a few miles of shore to 
offload the supply carriers. These supply carriers, LSVs, LCACs, etc., then must 
navigate through the mine fields as seen in Figure 4. This scenario would put a 
significant number of lives at risk and could easily lead to the loss of equipment, vessels, 
and personnel. Using the ASCC, the logistics vessel remains well outside of weapons 
range while offloading the ASCCs. The container would then travel along a defined route 
until it reaches its designated offload area. Once onshore, the ASCC could be unloaded 
manually, loaded on to cargo trucks for further transport, or reloaded with retrograde to 
return to the logistic ship off shore. Figure 5 shows the ASCCs route bypassing the mine 
field and supplying the troops on the beach with needed supplies. 



Figure 3. The ASCC Operating Environment (From: 14) 
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Figure 4. The Threat to Personnel Using Traditional Methods (From: 14) 



Figure 5. The ASCC Delivery System (From: 14) 


3. The AEPCO Design 

Advanced Engineering & Planning Corporation, Inc. (AEPCO) has proposed a 
Autonomous Sustainment Cargo Container that would solve the aforementioned 
problems. The ASCC aims to reduce manpower requirements for various littoral vessels 
as well as enhance personal safety by transporting cargo containers from the logistics 
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ship to shore using an autonomous system attached directly to the container. The typical 
ISO container can be seen in Figure 6 and comes in a 20 foot and 40 foot length. Initial 
designs have used to 20 foot ISO container to produce the vessel in Figure 7. The 
concept has a transit speed between 15 and 20 knots and enough range to transport goods 
from the logistics ship offshore, follow a programmed route past hazards, and arrive for 
pick up by troops ashore. AEPCO is currently involved in concept exploration in order to 
answer the fuel/power and speed/range questions. 



Figure 6. Typical ISO Container (From: 12) 


Ballasting Module 



Propulsion Ss'stcni 

Figure 7. AEPCO Proposed ASCC system (From: 13) 

E. THE PURPOSE OF THIS ANALYSIS 

The purpose of this analysis is to perform a proof of concept on the AEPCO 
ASCC design with an emphasis on seakeeping. The ASCC parameters will be input into 

a computer simulation/calculation program and the results will be compared with 
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indicated ship performance parameters. Based on the results, areas of concern or design 
shortfalls will be further scrutinized with an eye toward making design improvement 
suggestions. The simulated model will test the vessel in the areas of pitch, roll, vertical 
acceleration, lateral acceleration, mast submergence, propeller racing due to broaching, 
and bow slamming. These results will then be compared to probable operational areas to 
determine the operational limits to the particular operating environment. This will then 
allow for the calculation and display of the overall operational capability of this vessel. It 
is not the intention of this study to conduct a detailed structural analysis of the ASCC or 
create a finite element model of the seas and their interaction with the vessel structure. 
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II. SEAKEEPING ANALYSIS 


A. INITIAL SETUP 

I. The Proposed Layout 

The AEPCO ASCC is designed to use a standard 20 foot ISO container. It 
focuses on the 80% loading condition. This is both the most common and most 
economical way to deliver goods. This study will also use the 80% as the basis for 
analysis but will also add the 50, 60, and 70% loading conditions for comparison. 
Dimensions for the 20 foot ISO at 80% load can be seen in Table 1. This data was used 
to build the spreadsheet seen in Table 2. The ASCC KG calculation table is designed to 
solve for the various loading conditions based on a high, mid-box, or low loaded center 
of gravity. The distance from the keel to the center of gravity is called KG. The most 
common loading condition is the mid-box load where items inside the container are 
uniformly distributed and the center of gravity remains roughly in the center of the load. 
For ease of calculations, KG values were rounded to the nearest 10 for input into the 
analysis program. For the engine intake, a snorkel system has been designed to also act 
as a mast head to comply with Coast Guard Regulations. The mast head height is 
determined by the distance from which the mast head light can be seen. Coast Guard 
regulations require that partially submerged, unmanned vessels have a range of visibility 
of three miles.[8] This range converts to a required minimum mast head height of one 
meter. The one meter tall snorkel height is used for determining the mast submergence 
event occurrences. 
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Table 1. Dimensions of a 20 ft ISO Container at 80% Load (From: 13) 


1 20' ISO Under Loading Condition 

% Loaded 

80 

A [lbs] 

42316 

Length (L) [ft] 

19.875 

Beam (B) [ft] 

8 

Height (H) [ft] 

8 

A[LT] 

18.891 

▼ [ft^] 

661.188 

Draft (T) [ft] 

4.158 

Length-to- 

Beam Ratio 

(L/B) 

2.48 

Beam-to-Draft 

Ratio (B/T) 

1.92 

L' [ft] 

35.875 

L'/B 

4.48 
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Table 2. ASCC KG calculation 



Percent Load 

High Load 



0 

5 

10 

20 

30 

40 

50 

60 

70 

80 

Midbox Load 


weiqht 

8002.78 

10147.36 

12291.93 

16581.09 

20870.24 

25159.39 

29448.55 

33737.70 

38026.85 

42316.00 

Low Load 


weight (LT) 

3.57 

4.53 

5.49 

7.40 

9.32 

11.23 

13.15 

15.06 

16.98 

18.89 



V 

125.04 

158.55 

192.06 

259.08 

326.10 

393.12 

460.13 

527.15 

594.17 

661.19 



T 

0.79 

1.00 

1.21 

1.63 

2.05 

2.47 

2.89 

3.32 

3.74 

4.16 



KB 

0.39 

0.50 

0.60 

0.81 

1.03 

1.24 

1.45 

1.66 

1.87 

2.08 



BM 

54.25 

42.79 

35.32 

26.19 

20.80 

17.26 

14.74 

12.87 

11.42 

10.26 


KG (est) 

GM 












4.25 


50.40 

39.04 

31.68 

22.75 

17.58 

14.24 

11.94 

10.28 

9.04 

8.09 


4 


50.65 

39.29 

31.93 

23.00 

17.83 

14.49 

12.19 

10.53 

9.29 

8.34 


3.75 


50.90 

39.54 

32.18 

23.25 

18.08 

14.74 

12.44 

10.78 

9.54 

8.59 


3.5 


51.15 

39.79 

32.43 

23.50 

18.33 

14.99 

12.69 

11.03 

9.79 

8.84 


3.25 


51.40 

40.04 

32.68 

23.75 

18.58 

15.24 

12.94 

11.28 

10.04 

9.09 


3 


51.65 

40.29 

32.93 

24.00 

18.83 

15.49 

13.19 

11.53 

10.29 

9.34 


2.75 


51.90 

40.54 

33.18 

24.25 

19.08 

15.74 

13.44 

11.78 

10.54 

9.59 


2.5 


52.15 

40.79 

33.43 

24.50 

19.33 

15.99 

13.69 

12.03 

10.79 

9.84 


2.25 


52.40 

41.04 

33.68 

24.75 

19.58 

16.24 

13.94 

12.28 

11.04 

10.09 


2 


52.65 

41.29 

33.93 

25.00 

19.83 

16.49 

14.19 

12.53 

11.29 

10.34 


1.75 


52.90 

41.54 

34.18 

25.25 

20.08 

16.74 

14.44 

12.78 

11.54 

10.59 


1.5 


53.15 

41.79 

34.43 

25.50 

20.33 

16.99 

14.69 

13.03 

11.79 

10.84 


1.25 


53.40 

42.04 

34.68 

25.75 

20.58 

17.24 

14.94 

13.28 

12.04 

11.09 


1 


53.65 

42.29 

34.93 

26.00 

20.83 

17.49 

15.19 

13.53 

12.29 

11.34 



Metric 

Percent Load 


0 

5 

10 

20 

30 

40 

50 

60 

70 

80 

weiqht 

8002.78 

10147.36 

12291.93 

16581.09 

20870.24 

25159.39 

29448.55 

33737.70 

38026.85 

42316.00 

weight (LT) 

3.57 

4.53 

5.49 

7.40 

9.32 

11.23 

13.15 

15.06 

16.98 

18.89 

weight (kg) 

3630.06 

4602.84 

5575.62 

7521.18 

9466.74 

11412.30 

13357.86 

15303.42 

17248.98 

19194.54 

V(m") 

3.54 

4.49 

5.44 

7.33 

9.23 

11.13 

13.02 

14.92 

16.81 

18.71 

T(m) 

0.24 

0.30 

0.37 

0.50 

0.62 

0.75 

0.88 

1.01 

1.14 

1.27 

KB 

0.12 

0.15 

0.18 

0.25 

0.31 

0.38 

0.44 

0.50 

0.57 

0.63 

BM 

5.04 

3.98 

3.28 

2.43 

1.93 

1.60 

1.37 

1.20 

1.06 

0.95 

KG (est) 

GM 











1.2954 


53.35 

41.99 

34.63 

25.70 

20.53 

17.20 

14.90 

13.23 

11.99 

11.04 

1.2192 


53.43 

42.07 

34.71 

25.78 

20.61 

17.27 

14.97 

13.31 

12.07 

11.12 

1.143 


53.50 

42.14 

34.78 

25.86 

20.69 

17.35 

15.05 

13.38 

12.14 

11.20 

1.0668 


53.58 

42.22 

34.86 

25.93 

20.76 

17.43 

15.12 

13.46 

12.22 

11.27 

0.9906 


53.66 

42.30 

34.94 

26.01 

20.84 

17.50 

15.20 

13.54 

12.30 

11.35 

0.9144 


53.73 

42.37 

35.01 

26.09 

20.91 

17.58 

15.28 

13.61 

12.37 

11.43 

0.8382 


53.81 

42.45 

35.09 

26.16 

20.99 

17.66 

15.35 

13.69 

12.45 

11.50 

0.762 


53.88 

42.52 

35.16 

26.24 

21.07 

17.73 

15.43 

13.76 

12.52 

11.58 

0.6858 


53.96 

42.60 

35.24 

26.31 

21.14 

17.81 

15.50 

13.84 

12.60 

11.65 

0.6096 


54.04 

42.68 

35.32 

26.39 

21.22 

17.88 

15.58 

13.92 

12.68 

11.73 

0.5334 


54.11 

42.75 

35.39 

26.47 

21.30 

17.96 

15.66 

13.99 

12.75 

11.81 

0.4572 


54.19 

42.83 

35.47 

26.54 

21.37 

18.04 

15.73 

14.07 

12.83 

11.88 

0.381 


54.27 

42.90 

35.54 

26.62 

21.45 

18.11 

15.81 

14.15 

12.91 

11.96 

0.3048 


54.34 

42.98 

35.62 

26.69 

21.52 

18.19 

15.89 

14.22 

12.98 

12.03 


2. Modeling Waves and Sea States 

Vessel motions in regular waves were found using the strip theory of Salveson, 
Tuek, and Faltinsen as modified by Beek to inelude the surge degree of freedom.[3] The 
linear simultaneous equation that must be solved is 

,/=l-6 

^ +A.^) + ico^Bj^ +Cj^)g,^ =FI +F° where ^1-6 are surge, sway, heave, roll, 

1 = 1,6 

pitch, and yaw respectively; M is the mass matrix, A is the added mass matrix, B is the 
damping matrix, C is the hydrostatic restoring force matrix, is the Froude 
Krylovexciting force in the j* mode of motion, and F° is the diffraction exciting force in 
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the j*’’ mode of motion. For a port and starboard symmetrie ship the mass matrix is given 
as 

MO 0 0 Mz; 0 

0 M 0 -Mz, 0 Mx^ 

0 0 M 0 -Mx 0 

Mit = 

0 -Mz, 0 0 

Mz; 0 -My 0 y, o 

0 My, 0 0 4 

where M is the vessels total mass, x, y, and z are the eenters of mass in the longitudinal, 
transverse and vertieal directions, and I is the mass moment of inertia about the 
subscripted axis. In long crested random seas the response spectrum is determined by the 

linear relationship S^+{co,P) = \H{a),P)'^ Sj+ico) where Si+(ro) is the one-sided 

incident wave spectrum, H(ro,P) is the response amplitude operator for waves propagating 
in direction P, and So+(co,P) is the one-sided response spectrum due to waves traveling in 
direction p. A “cosine-squared” spreading function is used such that the area under the 
spectrum is the mean square value of the process. Several different wave spectra were 
used by the modeling program. The Pierson-Moskowitz spectrum as 

5^+( 0 ) = O.OOSlg^^ ^ where Hs is the significant wave height, ISSC 
spectrum as Sj+{a)) = i).\\Hla)'ll ^ where a)^=2nn^ and 

= mean period, ITTC spectrum based on zero crossings as 

Sj+{o)) = {HI !47r){co‘l / co ^) where co^ =lnI and 

= zero crossing period, ITTC spectrum based on modal (or peak) period as 

Sj + {(d) =)3\25HIo/^ /^ where co^ = 2;r /and 

= period of spectral peak, the Ochi 6-parameter spectrum as 

5, +(®) = 1/4X [{0.25(41. +1)<}1, 

J=l,2 

where = modal frequency of jth hump, 1^. = shape parameter for jth hump , and 
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li^j= significant wave height of jth hump , and 


JONS WAP spectrum as 


Sj +(®) = ag‘ 


CO 




\.15(o}^lmT) ^ /{2aX)) 


X=gxlUl 

f/jQ = wind speed at 10 meters 

X = fetch in meters 

, a = 0.076A 
where 

co„=2K3.5X^-^\glU,,) 


r = 3.3 

a = 0.07 for co < co^ 
<j = 0.09 for CO > co^^ 


Further refinement of the wave spectrum is done using the defined sea state 
significant wave height and modal wave period found in Table 3. The region specific 
percent probability values are used for later operability calculations. 


Table 3. Sea State Definition and Example Region Probabilities 


Sea 

State 

Significant Wave 
Height (m) 

Modal Wave 
Period (sec) 

Percent Probability 
(Region specific) 

1 

0.1 

6.3 

1 

2 

0.4 

6.5 

6 

3 

1.1 

7.5 

15 

4 

2.4 

8.8 

32 

5 

3.5 

9.7 

21 

6 

5.5 

12.4 

15 

7 

8.0 

15.0 

8 

8 

12.0 

16.5 

2 


3. Standards of Measurement for Naval Vessels 

The ASCC parametric limits were based on a combination of Naval standards 
shown in Figure 8, operational requirements, and mechanical limitations. Due to the 
unique nature of the ASCC, its lack of crew, radar, and sensitive electronics, motion 
extremes such as pitch, roll, slam, vertical acceleration and acceleration, and lateral 
velocity do not need to remain as fixed as the values listed in Figure 8. Pitch, roll, 
vertical and lateral accelerations, slam are generally associated with crew comfort. These 
motions contribute to seasic kn ess. For the ASCC, seasickness in the crew is not an issue 
and the parametrics can be freely modified to reflect the unique nature of this vessel. 
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These malleable limits were adjusted as follows: pitch limits were increased to 10 
degrees, roll went to 15 degrees, slam became 1.0 g’s, vertical acceleration was increased 
to 1.0 g’s, vertical velocity went to 3.0 m/s, and lateral accelerations became 0.5 g’s. For 
propeller racing events Figure 8 shows a limit of 25/100 which when calculated for a 
design trip of 40 NM leads to a value of 350 racing events per hour. A new category, 
mast submergence events, was added due to the damage that could occur should the 
engine intake pull a slug of water 10 cm or greater in length. For this study, the limit was 
set to the most conservative value, one event per trip. 
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Figure 8. 12 Prescribed Values for Naval Vessel Performance (From: 1) 

B. DESIGN ANALYSIS 
1. SHIPMO 

The program used to model the waves and wave responses is called 

SHIPMO.BM. The program predicts ship motions in six degrees of freedom (6-DOF) 

and five components of the shear and bending moment distributions. It is based on the 

strip theory approach of Salvesen, Tuck, and Faltinsen which was extended by Beck to 

include the surge degree of freedom. The program computes both regular and irregular 

waves and uses a cosine-squared spreading function for short crested seas. The program 

is based on the slender body theory, which assumes a long slender vessel with fine ends. 
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The version used for this analysis has been updated to inelude barge type forms with 
blunt ends whieh take into aeeount the end effeets. This is of partieular note for the box 
shaped ASCC. 9,196 runs were made using loads from 50 to 80 pereent in 10 pereent 
intervals, sea states from 1 to 8, and design speeds from 0 to 20 knots. 

a. Data Input 

Vessel speeifie data and run dynamies are input in a Fortran 77 format. 
The ships length, breadth, and design draft in meters are entered in at line 3. The vessels 
hull form is entered as a series of diserete points as seen in Figure 9. KG in meters 
measured from the waterline is entered as a difference with draft in line 8. The range of 
speeds in meters per second is entered into line 9. The range of headings in degrees for 
seas is in line 11. A sample input file can be found in Appendix B. Every change in 
draft, KG, and sea state becomes a new input fide, which leads to a large number of 
required runs to cover all likely combinations. 

68 (2FI0 4) {VAXIS<I,IS), ZAXIS(l,ISk l-f NT) 

Ther# ts on« d»te cord 6* for each input potnt Stert ot keel end work up 


Figure A-1 
input uffset Points 


The meximum i nput points 
(NT) l3 IS per etotlon 

Figure 9. Input Points Format for SHIPMO (From: 2) 

b. Data Extraction 

The process through which the SHIPMO program calculates the modeling 
results can be seen in the Computer Program Flow Chart in Figure 10 and 11. Each run 
in SHIPMO is output as a *.txt table and *.emf graphics file. These outputs are imported 
and collated through a previously created Matlab code shown in Appendix A. This code 
pulls together the sea state and load results and builds a graphic representation as 

17 




compared to the input ship limits as determined by ship standards in part A3 and mission 
requirements. The results eome out as a polar plot with ship heading toward the polar 
left, so waves striking the ships bow will be seen eoming from the 270 degree point on 
the polar plot. Speed rings are measured from the eenter being zero, the A ring at 5 
knots, the 14 ring at 10 knots, the 14 ring at 15 knots, and the outer ring at the max 
20knots. Measured results that exeeed the design parameters are noted in red, the various 
other shades note degrees to which the vessel is approaching those limits. 


Computer Program Flow Chart 



Figure 10. The SHIPMO Computer Program Flow Chart (From: 2) 
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STOP 


Figure 11. The SHIPMO Computer Program Flow Chart cont. (From: 2) 

C. RESULTS 

1. Initial Findings 

Preliminary review of the outputs showed a startling number of red areas in the 15 
to 20 knot speed range in head and fore quartering seas. Seleeted results ean be found in 
Appendix B. The output results for piteh, roll, slam, mast submergenee, propeller raeing, 
and lateral aeeeleration were review for potential problems in the design. The results 
were also put together to help determine what eategories were of more partieular 
importanee then others. Many of the results are based on the prineiples of ship motion 
involving erew eomfort and operational windows for erew laded vessels. Seasiek eausing 
motions like lateral and vertieal aeeelerations ean be eompletely ignored due to the laek 
of erew on the ASCC. Motions whieh eause seasickness but also have the potential to do 
damage, like slamming, can also be shown to be outside the scope of this study as the 
structural strength or placement of delicate components are not included. Propeller 
racing events although significant are less of a concern as the location of the propeller 
was not known and the data was estimated for a propeller 20 and 30 cm from the keel. 
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Additionally, Caterpillar diesels are extremely hardened against this type of marine 
environment abuse. After careful review, pitch, roll and mast submergence showed 
themselves as the most particularly significant results obtained. Mast submergence is 
especially important due to the possible damage that water intrusion can do to the 
Caterpillar diesel [15]. Both pitch and roll motions can be shown to be the root cause to 
both propeller racing and mast submergence events and are the driving factors to the 
limits to the ASCC’s operational capabilities. 

2. Results of Particular Importance 

a. Mast Submergence 

Based on documentation provided by Caterpillar Marine, the proposed 
supplier of the ASCCs prime mover, water intrusion of any amount into the intake would 
result in serious damage [15]. This damage is not relegated to simply the intake but 
could damage the manifold, turbo, exhaust and the block itself Additionally, a shortage 
of combustion air due to a water slug in the intake tube would lead to loss of combustion 
air and a stalling of the engine. Without crew onboard the ASCC would suffer this 
damage, stall, and be unable to continue its mission. This study took this particular event 
into account for analysis. A single mast submergence event on a 50 NM trip was 
considered unacceptable and can be seen in red in the following figures. Bow and fore 
quarter seas produce the worst results. Sea states above 3 in all configurations showed at 
least one event per trip at transit speeds over 15 knots. 



Figure 12. Mast Submergence at 80% Load, Sea State 1, Midbox Load. 
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Figure 13. 


Figure 14. 


Figure 15. 



Mast Submergence at 80% Load, Sea State 2, Midbox Load. 



Most Sub. (/trip); Loading ~ 80 %; KG ” 90 CTi; SS 3 Range ~ 1 NM; Maot ” 


270 


100 cm 

09 
0.8 
0.7 
0,6 
0.5 
0.4 
03 
0,2 
0,1 


Mast Submergence at 80% Load, Sea State 3, Midbox Load. 


Mast Sub- (/trip); Loading = 80 %; KG = 90 cn; SS 4 Range = 1 NM; Mast = 100 cm 



Mast Submergence at 80% Load, Sea State 4, Midbox Load. 
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Figure 16. 


Figure 17. 


Figure 18. 


Mast Sub. (/trip); Loading = 80 %; KG = 90 cn; SS 5 Range = 1 NM; Mast = 100 cm 

09 
0.8 
0.7 
0,6 
0.5 
0.4 
0.3 
0.2 
0.1 



Mast Submergence at 80% Load, Sea State 5, Midbox Load. 



Most Sub. (/trip); Loading ~ 80 %; KG ” 90 cm; SS 6 Range ~ 1 NM; Maot ” 


270 


100 cm 

09 
0.8 
0.7 
0,6 
0.5 
0.4 
03 
0.2 
0.1 


Mast Submergence at 80% Load, Sea State 6, Midbox Load. 


Mast Sub- (/trip); Loading = 80 %; KG = 90 cm; SS 7 Range = 1 NM; Mast = 100 cm 

09 
0.8 
0.7 
0.6 
0,5 
0 4 
0.3 
0.2 
0.1 



Mast Submergence at 80% Load, Sea State 7, Midbox Load. 
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Figure 19. 


Mast Submergence at 80% Load, Sea State 8, Midbox Load. 


b. Propeller Racing 

As previously noted, the propeller racing events are noteworthy simply to 
show the coincidence of pitch and roll to overall ships operability. As with mast 
submergence, bow and fore quartering seas produce the most events and these start as 
low as sea state 3. 



Figure 20. Propeller racing, 30 cm above keel at 80% Load, Sea State 1, Midbox 

Load. 
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Figure 21. 


Figure 22. 


Figure 23. 


Propeller Racing (/hr); Loading = SH %; KG = 90 cm; SS 2 ; Prop. = 30 cm 



Propeller racing, 30 cm above keel at 80% Load, Sea State 2, Midbox 

Load. 



30 cm 

I” 

- - 250 

• • 200 

■ ■ 150 

100 


Propeller Racing (/hr); Loading = 80 %; KG = 9C cm; SS 3 ; Prep. = 


270 



Propeller racing, 30 cm above keel at 80% Load, Sea State 3, Midbox 

Load. 


Propeller Racing (/hr); Loading = 80 %; KG = 9C cm; SS 4 ; Prep. = 30 cm 



Propeller racing, 30 cm above keel at 80% Load, Sea State 4, Midbox 

Load. 
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Figure 24. 


Figure 25. 


Figure 26. 


Propeller Racing (/hr); Loading = 80 %; KG = 90 cm; SS 5 ; Prop. = 30 cm 



Propeller racing, 30 cm above keel at 80% Load, Sea State 5, Midbox 

Load. 



Propeller racing, 30 cm above keel at 80% Load, Sea State 6, Midbox 

Load. 



Propeller racing, 30 cm above keel at 80% Load, Sea State 7, Midbox 

Load. 
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Prcpeller Racing (/hr); Loading = 80 %; KG = 90 cm; SS 8 ; Prop. = 30 cm 



Figure 27. Propeller racing, 30 cm above keel at 80% Load, Sea State 8, Midbox 

Load. 

By lowering the estimated propeller depth by 10 cm, the following figures 
are output. Although not readily obvious in this format, the change made had little effect 
on the likelihood of racing events. This will be addressed more thoroughly in a follow-on 
operability study. 



Figure 28. Propeller racing, 20 cm above keel at 80% Load, Sea State 1, Midbox 

Load. 
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Figure 29. 


Figure 30. 


Figure 31. 



Propeller racing, 20 cm above keel at 80% Load, Sea State 2, Midbox 

Load. 


Propeller Racing (/hr); Loading = 80 %; KG = 9C cm; SS 3 ; Prep. = 20 cm 



Propeller racing, 20 cm above keel at 80% Load, Sea State 3, Midbox 

Load. 

Propeller Racing (/hr); Loading = 83 %; KG = 90 cm; SS 4 ; Prop. = 20 cm 



Propeller racing, 20 cm above keel at 80% Load, Sea State 4, Midbox 

Load. 
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Figure 32. 


Propeller Racing (/hr); Loading = 83 %; KG = 90 cm; SS 5 ; Prop. = 20 cm 



Propeller racing, 20 cm above keel at 80% Load, Sea State 5, Midbox 

Load. 
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Figure 33. 


Figure 34. 


Prcpcilcr Racing (/hr); Loading - 80 %; KG - 90 cm; SS 6 ; Prop. “ 20 cm 

300 

250 

2UU 

150 

100 

50 

270 I 

Propeller racing, 20 cm above keel at 80% Load, Sea State 6, Midbox 

Load. 




270 


Propeller Racing (/ht); Loading = 80 %; KG = 90 :m; SS 7 ; Prop. 


250 

230 

150 

130 

53 


Propeller racing, 20 cm above keel at 80% Load, Sea State 7, Midbox 

Load. 
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Figure 35. Propeller racing, 20 cm above keel at 80% Load, Sea State 8, Midbox 

Load. 

c. Pitch and Roll 

The results show a distinet interaetion between piteh and roll and mast 
submergenee events and propeller raeing events. Note in the following figures that both 
pitch and roll degrees go in extremis with waves on the bow or the fore quarters. It is 
also noteworthy that where the red areas of the pitch and roll motions overlap, both mast 
submergence and propeller racing go red. 

Pitch (deg) Loading = 80 %; KG = 90 cm; SS 1 

0.9 
U.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

270 I _I 



Figure 36. Pitch at 80% Load, Sea State 1, Midbox Load. 
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Figure 37. 


Figure 38. 


Figure 39. 


Pitch (deg) Loading = 80 %; KG = 90 cm; SS 2 



Pitch at 80% Load, Sea State 2, Midbox Load. 


Pitch (deg); Loading = 80 %; KG = 90 cm; SS 3 



Pitch at 80% Load, Sea State 3, Midbox Load. 


Pitch (deg) Loading = 80 %: KG = 90 cm; SS 4 



Pitch at 80% Load, Sea State 4, Midbox Load. 
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Figure 40. 


Figure 41. 


Figure 42. 


Pitch (deg); Loading = 80 %; KG = 90 cm; SS 5 



Pitch at 80% Load, Sea State 5, Midbox Load. 


Pitch (deg); Loading = 80 %; KG =90 cm; SS 6 



Pitch at 80% Load, Sea State 6, Midbox Load. 


Pitch (deg); Loading = 80 %; KG = 93 cm; SS 7 



Pitch at 80% Load, Sea State 7, Midbox Load. 
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Figure 43. 


Figure 44. 


Figure 45. 


Pitch (deg); Loading = 80 %; KG = 90 cm; SS 8 



Pitch at 80% Load, Sea State 8, Midbox Load. 


Roll (deg); Loading = 80 %. KG =90 cm; SS 1 



Roll at 80% Load, Sea State 1, Midbox Load. 


Roll (deg); Loading = 80 %\ KG =90 cm; SS 2 



Roll at 80% Load, Sea State 2, Midbox Load. 
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Figure 46. 


Figure 47. 


Figure 48. 


Roll (deg); Loading = 80 %; KG = 90 cm; S3 3 



Roll at 80% Load, Sea State 3, Midbox Load. 


Roll (deg); Loading = 80 %; KG = 90 cm; SS 4 



Roll at 80% Load, Sea State 4, Midbox Load. 


Roll (deg); Loading = 80 %\ KG =90 cm; SS 5 



Roll at 80% Load, Sea State 5, Midbox Load. 
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Figure 49. 


Figure 50. 


Figure 51. 


Roll (deg); Loading = 80 %; KG = 90 cm; S3 6 



Roll at 80% Load, Sea State 6, Midbox Load. 



Roll (deg); Loading = 80 %; KG = 90 cm; SS 7 


270 



Roll at 80% Load, Sea State 7, Midbox Load. 


Roll (deg); Loading = 80 %\ KG =90 cm; SS 8 



270 


Roll at 80% Load, Sea State 8, Midbox Load. 
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D. CONCLUSSIONS 


1. Design/Mission Shortcomings 

While the goal of the ASCC is to operate in sea states above the current systems, 
Figures 52 and 53 indicate that this may not be entirely possible. By putting in a regional 
probability of occurrence for the different sea states and overlaying them with the 
operability of the vessel with all of the previously mentioned parametrics, it can be seen 
that only in sea states 1 and 2 is there a low occurrence of operational interruptions. At 
sea state 3 the ASCC would only be able to operate less than 65% of the time. At sea 
state 4, the set goal, operability drops below 45%. It must be noted that these figures are 
for all speeds and headings. It would be possible to increase the operability of the ASCC 
by mandating wave heading and slower speeds, but under the current design sea state 4 is 
simply too damaging to be able to successfully operate. Also note that very little 
difference is made between the results of the two assumed propeller height values. 


Loading = £0 %; KG = 90 cm Range = 1 NM; Mast = 100 cm; Prop. = 30 cm 



1 2 3 4 5 6 7 8 

Sea State 


Figure 52. Overall Operability Index for 80% Midbox Load and 30 cm propeller 

height. 
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Loading = £0 %; KG= 90 cm Range = 1 NM; Mast = 100 cm; Prop. = 20 cm 



1 2 3 4 5 6 7 8 

Sea State 


Figure 53. Overall Operability Index for 80% Midbox Load and 20 cm propeller 

height. 

2. Improvement Possibilities 

Although the predicted response is deemed satisfactory for the intended purpose 
of the ASCC, it seakeeping performance could be improved as follows: (1) Increased 
mast height in order to reduce mast submergence events. Shown below are mast 
submergence events for mast height configurations from 1 to 2 meters, 50% loaded 
(return load), a KG of 50 cm (upper midbox), and at the design goal of sea state 4. As 
can bee seen, raising the mast head height and engine intake significantly reduces the 
number of submergence events down to very rare for the 2 m case. This, however, might 
increase the center of gravity of the system and would need further testing. (2) Increased 
roll damping through rudder control. This, however, might increase the added resistance 
due to rudder action and, therefore, reduce the available payload. (3) Ballasted outriggers 
in order to decrease roll period. This, however, might increase resistance and complicate 
launch and recovery, and limit the number of containers that can be carried by the 
logistics ship. 
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Figure 54. 


Figure 55. 


Figure 56. 


Mast Sub. (/trip); Loading = 50 %; KG = 50 cm; SS 4 ; Range = 4 NM; Mast = 100 cm 

D 9 

a.8 

0.7 
0.6 
0,5 
0 4 
0.3 
0.2 
0.1 



Mast submergence at the initial mast height of 1 m. 



Mast Sub. (/tnp); Loading = 50 %; KG = 50 cm; SS 4 , Range = 4 NM; Mas1 = 
90 in 


120 cm 

0.9 
0.8 


0.3 

0.2 



Mast submergence with a mast height of 1.2 m. 



Mast submergence with a mast height of 1.4 m. 
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Figure 57. 


Figure 58. 


Figure 59. 


Mast Sub. (/trip); Loading = 50 %; KG = 50 cm; SS 4 ; Range = 4 NM; Mast = 160 cm 

D.9 
□ .8 
0,7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



Mast submergence with a mast height of 1.6 m. 


Mast Sub. i7trip); Loading = 50 %; <G = 50 cm; SS 4 ; Range = 4 NM; Mast = 180 cm 



Mast submergence with a mast height of 1.8 m. 


Mast Sub. (/trip): Loading = 53 %; KG = 50 cm; SS 4 ; Range = 4 NVI; Mas! = 200 cm 



Mast submergence with a mast height of 2 m. 
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3. The Next Step 

More detailed seakeeping ealeulations eould be performed by taking into 
eonsideration the exact shape of the bow and stem units, along with experimental 
verification. Detailed graphical modeling would also be helpful in further determining 
root causes to the mast submergence and propeller racing events. Agencies such as the 
David Taylor Institute and Naval Surface Warfare Center Carderock could be engaged to 
conduct both digital FEM and CDF modeling and physical tank model experiments to 
further refine these results. 
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APPENDIX A. MATLAB COLLATION LUNCTION 


function fcn_polar_speed_loop(Loading_Cond,KG,z_prop_keel,z_mast_top,range_NM) 
% ASCC Container Data Set 

% Horizontal/Vertical motions - Polar plots 
% Speed in knots - Heading (180 head, 0 stern) 

% Contour plots (heading/speed) 

% Two parameter Bretchneider spectrum - Long crested seas 
% Dimensional version (Metric units) 


% General Loop - meant for overnight runs 
warning off 
% GENERAL DATA 


lambda_min 

lambda_max 

delta_lambda 

rho 

zeta 

L 

g 


= 5.000; 

= 505.000; 
= 10.060; 
=1024.000; 
1 . 000 ; 
6.058; 
9.810; 


% Sea state information 


% Min wave length (m) 

% Max wave length (m) 

% Wave length increment (m) 

% Water density 
% Regular wave height 

% Reference length for nondimensionalization 
% Gravitational constant 


sig_wave_height 

modal_period_sec 

percent_probab 


[0.1 0.4 1 
[6.3 6.5 7 
[1 6 15 32 


12.43.5 5.5 8.0 12.0]; 
5 8.8 9.7 12.4 15.0 16.5]; 
21 15 8 2]; 


% Set limits for operations (PNA Vol. Ill) 

% All values in RMS single amplitudes 
% at Point-of-Interest (POI) where applicable 


c_p_roll_max = 

15; 

c_p_pitch_max = 

o 
;—1 

c_p_ver_acc_max = 

O 

;—1 

c_p_lat_acc_max = 

0.5; 

c_p_bow_acc_max = 

o 

;—1 

c_p_bow_vel_max = 

3.0; 

c_p_prop_max = 

350; 

c_p_mast_trip_max 

= 


% Roll (deg) 

% Pitch (deg) 

% Abs . Vert. Accel, at amidships (g's) 
% Abs. Lat . Accel, at amidships (g's) 

% Slam Accel, at bow (g's) 

% Relative vert, velocity at bow (m/s) 
% Propeller racing events (per hour) 

% Mast submergence events (per trip) 


c_p_roll = 
c_p_pitch = 
c_p_ver_acc = 
c_p_lat_acc = 
c_p_bow_acc = 
c_p_b o w_V el = 
c_p_prop = 
c_p_mast_trip 


[0:0.5:c_p_roll_max]; 

[0:1:c_p_pitch_max] ; 

[0:0.05:c_p_ver_acc_max]; 
[0:0.02:c_p_lat_acc_max]; 
[0:0.05:c_p_bow_acc_max]; 
[0:0.10:c_p_bow_vel_max]; 
[0:10:c_p_prop_max]; 

= [0:0.1:c_p_mast_trip_max]; 


c_p_roll_max_string = 
c_p_pitch_max_string = 
c_p_ve r_ac c_max_s t ring = 
c_p_lat_acc_max_string = 
c_p_b o w_a c c_m a x_ string = 
c_p_b o w_V el_max_string = 
c_p_prop_max_string = 
c_p_mast_trip_max_string 


num2str(c_p_roll_max); 
num2str(c_p_pitch_max); 
num2str(c_p_ver_acc_max); 
num2str(c_p_lat_acc_max); 
num2str(c_p_bow_acc_max); 
num2str(c_p_bow_vel_max); 
num2str(c_p_prop_max); 

= num2str(c_p_mast_trip_max); 


% Speed range (m/sec) 


V_min = 0; 

V_max = 10; 
delta_V = 1; 

Speed_Range_string= 'All Speeds ' ; 
Loading_Cond_string = num2str(Loading_Cond); 
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if Loading_Cond == 80 
draft = 1.27; 

end 

if Loading_Cond == 70 
draft = 1.14; 

end 

if Loading_Cond == 60 
draft = 1.01; 

end 

if Loading_Cond == 50 
draft = 0.88; 

end 

KG_string = num2str(KG); 

% Point for motions calculation 

z_point = 0.00; 
x_point = 0.00; 
y_point = 0.00; 

z_prop_keel_string = num2str(z_prop_keel); 
z_prop = z_prop_keel/100 - draft; 
z_prop_string = num2str(z_prop); 
x_prop = -L/2; 
y_prop = 0.00; 

z_mast = z_mast_top/100 + 2.5 - draft; 
z_mast_top_string = num2str(z_mast_top); 
z_mast_string = num2str(z_raast); 
x_mast = -L/2; 
y_mast = 0.00; 

range_string = num2str(range_NM); 


% Loop over sea state 

i_Sea_State=0; 

for Sea_State=l:1:8, 

i_Sea_State=i_Sea_State+l; 

[Loading_Cond KG z_prop_keel z_mast_top range_NM Sea_State] 
S_S_string = num2str(Sea_State); 

HS=sig_wave_height(Sea_State); 

HS_string=num2str(HS); 

T_m=modal_period_sec(Sea_State); 

T_m_string = num2str(T_m); 
omega_m=2*pi/T_m; 


% Set up file reading format. 


trigg = 30; 
flloc = 25; 
f21oc = 26; 
f31oc = 27; 
f41oc = 28; 
f51oc = 29; 
f61oc = 30; 


lambda 

lambdaoverL 

wavenumber 

omega 

wavenumber 
omega 
period 
filesize 
lambda_size 


lambda_min:delta_lambda:lambda_max; 
lambda/L; 

2.0*pi./lambda; 
sqrt(wavenumber*g); 
wavenumber'; 
omega'; 

2.0*pi./omega; 
size(lambda); 
trigg*filesize(2); 


beta_incr = 10; 


D Vector of wavelengths 

% Wave number 
% Wave frequency 
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iV=0; 

for V = V_min:delta_V:V_max, 

iV = iV+1; 

V_string = num2str(V); 

V_calc = V; % m/sec 

V_knots = V/0.5144; 

tirae_hrs = range_NM/(V_knots+0.01); 

ibeta = 0; 

for beta = 0:beta_incr:360, % Loop on sea direction 

if beta > 180 

beta= 360-beta; 
end 

beta_string = num2str(beta); 
beta_calc = beta*pi/180; 
ibeta = ibeta+1; 

% Load data file V_beta.txt 


load_filename = 

strcat (' results\ ' ,Loading_Cond_string,,KG_string,,V_string, 
filename = load(load_filename); 


omegae = omega-wavenumber*V_calc*cos(beta_calc); 

periode = 2.0*pi./omegae; 


% Set mass matrix elements 


Mil 
M12 
Ml 3 
M14 
M15 
Ml 6 
M21 
M2 2 
M2 3 
M2 4 
M2 5 
M2 6 
M31 
M32 
M33 
M34 
M35 
M36 
M41 
M42 
M43 
M44 
M4 5 
M4 6 
M51 
M52 
M53 
M54 
M55 
M56 
M61 
M62 
M63 
M64 
M65 
M66 


filename (1: trigg: lambda, 
filename (1: trigg: lambda, 
filename (1: trigg: lambda, 
filename (1: trigg: lambda, 
filename (1: trigg: lambda, 
filename (1: trigg: lambda, 
filename (2 : trigg: lambda, 
filename (2 : trigg: lambda, 
filename (2 : trigg: lambda, 
filename (2 : trigg: lambda, 
filename (2 : trigg: lambda, 
filename (2 : trigg: lambda, 
filename (3 : trigg: lambda, 
filename (3 : trigg: lambda, 
filename (3 : trigg: lambda, 
filename (3 : trigg: lambda, 
filename (3 : trigg: lambda, 
filename (3 : trigg: lambda, 
filename (4 : trigg: lambda, 
filename (4 : trigg: lambda, 
filename (4 : trigg: lambda, 
filename (4 : trigg: lambda, 
filename (4 : trigg: lambda, 
filename (4 : trigg: lambda, 
filename (5 : trigg: lambda, 
filename (5 : trigg: lambda, 
filename (5 : trigg: lambda, 
filename (5 : trigg: lambda, 
filename (5 : trigg: lambda, 
filename (5 : trigg: lambda, 
filename (6 : trigg: lambda, 
filename (6 : trigg: lambda, 
filename (6 : trigg: lambda, 
filename (6 : trigg: lambda, 
filename (6 : trigg: lambda, 
filename (6 : trigg: lambda. 


[Mil 

M12 

M13 

M14 

M15 

Ml 6 

M21 

M2 2 

M2 3 

M2 4 

M2 5 

M2 6 

M31 

M32 

M33 

M34 

M35 

M36 

M41 

M42 

M4 3 

M4 4 

M4 5 

M4 6 

M51 

M52 

M53 

M54 

M55 

M5 6 

M61 

M62 

M63 

M64 

M65 

M66 


s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 

s 


ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize, 

ize. 


1 ) ; 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 


_ ' ,beta_string, '. txt ' ) 

Frequency of encounter 
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% Added mass terms 

All = filename{6+1:t 
A12 = filename{6+1:t 
A13 = filename{6+1:t 
A14 = filename{6+1:t 
A15 = filename{6+1:t 
A16 = filename{6+1:t 
A21 = filename{6+2:t 
A22 = filename{6+2:t 
A23 = filename{6+2:t 
A24 = filename{6+2:t 
A25 = filename{6+2:t 
A26 = filename{6+2:t 
A31 = filename{6+3:t 
A32 = filename{6+3:t 
A33 = filename{6+3:t 
A34 = filename{6+3:t 
A35 = filename{6+3:t 
A36 = filename{6+3:t 
A41 = filename{6+4:t 
A42 = filename{6+4:t 
A43 = filename{6+4:t 
A44 = filename{6+4:t 
A45 = filename{6+4:t 
A46 = filename{6+4:t 
A51 = filename{6+5:t 
A52 = filename{6+5:t 
A53 = filename(6+5:t 
A54 = filename{6+5:t 
A55 = filename(6+5:t 
A56 = filename(6+5:t 
A61 = filename{6+6:t 
A62 = filename{6+6:t 
A63 = filename{6+6:t 
A64 = filename{6+6:t 
A65 = filename(6+6:t 
A66 = filename(6+6:t 

A = [All A12 A13 A14 
A21 A22 A23 A24 
A31 A32 A33 A34 
A41 A42 A43 A44 
A51 A52 A53 A54 
A61 A62 A63 A64 


rigg: larrLbda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 
rigg:lambda_size, 

A15 A16;... 

A2 5 A2 6; . . . 

A35 A36; . . . 

A45 A46;... 

A55 A56;... 

A65 A66]; 


% Damping terms 


Bll 

= 

filename 

(12 + 1 

trigg 

B12 

= 

filename 

(12 + 1 

trigg 

B13 

= 

filename 

(12 + 1 

trigg 

B14 

= 

filename 

(12 + 1 

trigg 

B15 

= 

filename 

(12 + 1 

trigg 

B16 

= 

filename 

(12 + 1 

trigg 

B21 

= 

filename 

(12 + 2 

trigg 

B22 

= 

filename 

(12 + 2 

trigg 

B23 

= 

filename 

(12 + 2 

trigg 

B24 

= 

filename 

(12 + 2 

trigg 

B25 

= 

filename 

(12 + 2 

trigg 

B2 6 

= 

filename 

(12 + 2 

trigg 

B31 

= 

filename 

(12 + 3 

trigg 

B32 

= 

filename 

(12 + 3 

trigg 

B33 

= 

filename 

(12 + 3 

trigg 

B34 

= 

filename 

(12 + 3 

trigg 

B35 

= 

filename 

(12 + 3 

trigg 

B36 

= 

filename 

(12 + 3 

trigg 

B41 

= 

filename 

(12 + 4 

trigg 

B42 

= 

filename 

(12 + 4 

trigg 

B43 

= 

filename 

(12 + 4 

trigg 

B44 

= 

filename 

(12 + 4 

trigg 

B45 

= 

filename 

(12 + 4 

trigg 


lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 

lambda_size 


1 ) ; 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 
1 ); 
2 ) ; 

3) ; 

4) ; 

5) ; 

6 ) ; 


, 1 ); 
,2) ; 
,3) ; 
,4) ; 
,5) ; 
, 6 ); 
, 1 ) ; 
, 2 ) ; 
,3) ; 
,4) ; 
,5) ; 
, 6 ); 
, 1 ) ; 
, 2 ) ; 
, 3 ) ; 
,4) ; 
,5) ; 
, 6 ); 
, 1 ) ; 
, 2 ) ; 
, 3 ) ; 
,4) ; 
,5) ; 
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B46 = filename (12 + 4 : trigg: larrLbda_size, 6) ; 
B51 = filename (12 + 5 : trigg: larrLbda_size, 1) ; 
B52 = filename (12 + 5 : trigg: larrLbda_size, 2) ; 
B53 = filename(12+5:trigg:lambda_size,3); 
B54 = filename(12+5:trigg:lambda_size,4); 
B55 = filename(12+5:trigg:lambda_size,5); 
B56 = filename(12+5:trigg:lambda_size,6); 
B61 = filename(12+6:trigg:lambda_size,1); 
B62 = filename(12+6:trigg:lambda_size,2); 
B63 = filename(12+6:trigg:lambda_size,3); 
B64 = filename(12+6:trigg:lambda_size,4); 
B65 = filename(12+6:trigg:lambda_size,5); 
B66 = filename(12+6:trigg:lambda_size,6); 


[Bll 

B12 

B13 

B14 

B15 

B16;.. 

B21 

B22 

B23 

B24 

B25 

B26; . . 

B31 

B32 

B33 

B34 

B35 

B36; . . 

B41 

B42 

B43 

B44 

B45 

B46; . . 

B51 

B52 

B53 

B54 

B55 

B56; . . 

B61 

B62 

B63 

B64 

B65 

B66] ; 


% Hydrostatic terms 


Cll 

C12 

C13 

C14 

C15 

C16 

C21 

C22 

C23 

C24 

C25 

C26 

C31 

C32 

C33 

C34 

C35 

C36 

C41 

C42 

C43 

C44 

C45 

C46 

C51 

C52 

C53 

C54 

C55 

C56 

C61 

C62 

C63 

C64 

C65 

C66 


filename(18+l:trigg 
filename(18+1:trigg 
filename(18+1:trigg 
filename(18+1:trigg 
filename(18+1:trigg 
filename(18+l:trigg 
filename(18+2:trigg 
filename(18+2:trigg 
filename(18+2:trigg 
filename(18+2:trigg 
filename(18+2:trigg 
filename(18+2:trigg 
filename(18+3:trigg 
filename(18+3:trigg 
filename(18+3:trigg 
filename(18+3:trigg 
filename(18+3:trigg 
filename(18+3:trigg 
filename(18+4:trigg 
filename(18+4:trigg 
filename(18+4:trigg 
filename(18+4:trigg 
filename(18+4:trigg 
filename(18+4:trigg 
filename(18+5:trigg 
filename(18+5:trigg 
filename(18+5:trigg 
filename(18+5:trigg 
filename(18+5:trigg 
filename(18+5:trigg 
filename(18+6:trigg 
filename(18+6:trigg 
filename(18+6:trigg 
filename(18+6:trigg 
filename(18+6:trigg 
filename(18+6:trigg 


lambda_size,1); 
lambda_size,2); 
lambda_size,3); 
lambda_size,4); 
lambda_size,5); 
lambda_size,6); 
lambda_size,1); 
lambda_size,2); 
lambda_size,3); 
lambda_size,4); 
lambda_size,5); 
lambda_size,6); 
lambda_size,1); 
lambda_size,2); 
lambda_size,3); 
lambda_size,4); 
lambda_size,5); 
lambda_size,6); 
lambda_size,1); 
lambda_size,2); 
lambda_size,3); 
lambda_size,4); 
lambda_size,5); 
lambda_size,6); 
lambda_size,1); 
lambda_size,2); 
lambda_size,3); 
lambda_size,4); 
lambda_size,5); 
lambda_size,6); 
lambda_size,1); 
lambda_size,2); 
lambda_size,3); 
lambda_size,4); 
lambda_size,5); 
lambda_size,6); 


C = [Cll C12 C13 C14 C15 C16;. 
C21 C22 C23 C24 C25 C26;. 
C31 C32 C33 C34 C35 C36;. 
C41 C42 C43 C44 C45 C46;. 
C51 C52 C53 C54 C55 C56;. 
C61 C62 C63 C64 C65 C66]; 


% Total exciting forces 


Fl_amp = filename(fHoc:trigg:lambda_size,5); 
F2_amp = filename(f21oc:trigg:lambda_size,5); 
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ASlbar. 

Allbar. 
ASlbar. 

Allbar. 
A5 Ibar. 


A42bar. 

A22bar. 
A62bar. 

A22bar. 
A62bar. 


A31bar. 

Allbar. 
ASlbar. 

Allbar. 
AS Ibar. 


F 3_amp 
F4_amp 
FS_amp 
F 6_amp 


filename(f31oc:trigg:lambda_size,S); 
filename(f41oc:trigg:lambda_size,S); 
filename(fSloe:trigg:lambda_size,S); 
filename(f61oc:trigg:lambda_size,S); 


Fl_pha 

F2_pha 

F3_pha 

F4_pha 

FS_pha 

F6_pha 


filename(fHoc:trigg:lambda_size, 6) ; 
filename(f21oc:trigg:lambda_size,6); 
filename(f31oc:trigg:lambda_size,6); 
filename(f41oc:trigg:lambda_size,6); 
filename(fSloe:trigg:lambda_size,6); 
filename(f61oc:trigg:lambda_size,6); 


FI = Fl_amp.*exp(i*Fl_pha.*pi/180.0); 
F2 = F2_amp.*exp{i*F2_pha.*pi/180.0); 
F3 = F3_amp.*exp{i*F3_pha.*pi/180.0); 
F4 = F4_amp.*exp{i*F4_pha.*pi/180.0); 
FS = FS_amp.*exp{i*FS_pha.*pi/180.0); 
F6 = F6_amp.*exp{i*F6_pha.*pi/180.0); 

F = [F1;F2;F3;F4;FS;F6]; 

% Calculate motions 


Allbar 

= 

- {omegae.^2) .* 

(Mil 

+ 

All) 

+ 

i*omegae.*B11 

+ 

Cll 

A13bar 

= 

-{omegae .^ 2 ) .* 

(M13 

+ 

A13) 

+ 

i*omegae.*B13 

+ 

C13 

AlSbar 

= 

-{omegae .^ 2 ) . * 

(M15 

+ 

A15) 

+ 

i*omegae.*B1S 

+ 

CIS 

A22bar 

= 

- {omegae .'^ 2 ) . * 

(M22 

+ 

A22) 

+ 

i*omegae.*B22 

+ 

C22 

A24bar 

= 

-{omegae .^ 2 ) .* 

(M2 4 

+ 

A24) 

+ 

i*omegae.*B24 

+ 

C24 

A2 6bar 

= 

-{omegae .^ 2 ) .* 

(M2 6 

+ 

A2 6) 

+ 

i*omegae.*B2 6 

+ 

C2 6 

A31bar 

= 

-{omegae .^ 2 ) .* 

(M31 

+ 

A31) 

+ 

i*omegae.*B31 

+ 

C31 

A33bar 

= 

-{omegae .^ 2 ) .* 

(M33 

+ 

A33) 

+ 

i*omegae.*B33 

+ 

C33 

A3Sbar 

= 

-{omegae .^ 2 ) .* 

(M35 

+ 

A35) 

+ 

i*omegae.*B3S 

+ 

C35 

A42bar 

= 

-{omegae .^ 2 ) . * 

(M42 

+ 

A42) 

+ 

i*omegae.*B42 

+ 

C42 

A44bar 

= 

-{omegae .^ 2 ) .* 

(M44 

+ 

A44) 

+ 

i*omegae.*B44 

+ 

C44 

A4 6bar 

= 

-{omegae .^ 2 ) . * 

(M4 6 

+ 

A4 6) 

+ 

i*omegae.*B4 6 

+ 

C46 

ASIbar 

= 

-{omegae .^ 2 ) . * 

(M51 

+ 

A51) 

+ 

i*omegae.*BS1 

+ 

C51 

AS3bar 

= 

-{omegae .^ 2 ) .* 

(M53 

+ 

A53) 

+ 

i*omegae.*BS3 

+ 

C53 

ASSbar 

= 

-{omegae .^ 2 ) .* 

(M55 

+ 

A55) 

+ 

i*omegae.*BSS 

+ 

CSS 

A62bar 

= 

-{omegae .^ 2 ) .* 

(M62 

+ 

A62) 

+ 

i*omegae.*B62 

+ 

C62 

A64bar 

= 

-{omegae .^ 2 ) .* 

(M64 

+ 

A64) 

+ 

i*omegae.*B64 

+ 

C64 

A66bar 

= 

-{omegae.^2).* 

(M66 

+ 

A66) 

+ 

i*omegae.*B66 

+ 

C66 


etal_c = (A33bar.*ASSbar - A3Sbar.*AS3bar)./... 

{A1Ibar.*A33bar.*ASSbar - A1Ibar.*A3Sbar.*AS3bar - A31bar.*A13bar.*ASSbar + 
*AlSbar.*AS3bar + ASIbar.*A13bar.*A3Sbar - ASlbar.*AlSbar.*A33bar).*F1 - ... 

{A13bar.*ASSbar - AlSbar.*AS3bar)./{Allbar.*A33bar.*ASSbar - 
*A3Sbar.*AS3bar - A31bar.*A13bar.*ASSbar + A31bar.*AlSbar.*AS3bar + 

*A13bar.*A3Sbar - ASlbar.*AlSbar.*A33bar).*F3 + ... 

(A13bar.*A3Sbar - AlSbar.*A33bar)./(Allbar.*A33bar.*ASSbar - 
*A3Sbar.*AS3bar - A31bar.*A13bar.*ASSbar + A31bar.*AlSbar.*AS3bar + 

*A13bar.*A3Sbar - ASIbar.*AlSbar.*A33bar).*FS; 
eta2_c = (A44bar.*A66bar - A46bar.*A64bar)./... 

(A22bar.*A44bar.*A66bar - A22bar.*A46bar.*A64bar - A42bar.*A24bar.*A66bar + 
*A26bar.*A64bar + A62bar.*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F2 - ... 

{A24bar.*A66bar - A26bar.*A64bar)./(A22bar.*A44bar.*A66bar - 
*A46bar.*A64bar - A42bar.*A24bar.*A66bar + A42bar.*A26bar.*A64bar + 

*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F4 + ... 

(A24bar.*A46bar - A26bar.*A44bar)./(A22bar.*A44bar.*A66bar - 
*A46bar.*A64bar - A42bar.*A24bar.*A66bar + A42bar.*A26bar.*A64bar + 

*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F6; 
eta3_c = -(A31bar.*ASSbar - A3Sbar.*ASIbar)./... 

(A1Ibar.*A33bar.*ASSbar - A1Ibar.*A3Sbar.*AS3bar - A31bar.*A13bar.*ASSbar + 
*AlSbar.*AS3bar + ASlbar.*A13bar.*A3Sbar - ASlbar.*AlSbar.*A33bar).*F1 + ... 

{Allbar.*ASSbar - AlSbar.*ASlbar)./{Allbar.*A33bar.*ASSbar - 
*A3Sbar.*AS3bar - A31bar.*A13bar.*ASSbar + A31bar.*AlSbar.*AS3bar + 

*A13bar.*A3Sbar - ASlbar.*AlSbar.*A33bar).*F3 - ... 

{A1Ibar.*A3Sbar - AlSbar.*A31bar)./{A1Ibar.*A33bar.*ASSbar - 
*A3Sbar.*AS3bar - A31bar.*A13bar.*ASSbar + A31bar.*AlSbar.*AS3bar + 

*A13bar.*A3Sbar - ASIbar.*AlSbar.*A33bar).*FS; 
eta4_c = -{A42bar.*A66bar - A46bar.*A62bar)./... 
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A42bar. 


A22bar. 
A62bar. 

A22bar. 
A62bar. 


ASlbar. 

Allbar. 
ASlbar. 

Allbar. 
A5 Ibar. 


A42bar. 

A22bar. 
A62bar. 

A22bar. 
A62bar. 


(A22bar.*A44bar.*A66bar - A22bar.*A46bar.*A64bar - A42bar.*A24bar.*A66bar + 
*A26bar.*A64bar + A62bar.*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F2 + ... 

{A22bar.*A66bar - A26bar.*A62bar)./{A22bar.*A44bar.*A66bar - 
*A46bar.*A64bar - A42bar.*A24bar.*A66bar + A42bar.*A26bar.*A64bar + 

*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F4 - ... 

(A22bar.*A46bar - A26bar.*A42bar)./{A22bar.*A44bar.*A66bar - 
*A46bar.*A64bar - A42bar.*A24bar.*A66bar + A42bar.*A26bar.*A64bar + 

*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F6; 
eta5_c = {ASlbar.*A53bar - ASSbar.*A5Ibar)./... 

(A1Ibar.*A33bar.*A55bar - A1Ibar.*A35bar.*A53bar - ASlbar.*A13bar.*A55bar + 
*A15bar.*A53bar + A5Ibar.*A13bar.*A35bar - A5Ibar.*A15bar.*A33bar).*F1 - ... 

(Allbar.*A53bar - AlSbar.*A51bar)./(Allbar.*A33bar.*A55bar - 
*A35bar.*A53bar - ASlbar.*A13bar.*A55bar + ASlbar.*A15bar.*A53bar + 

*A13bar.*A35bar - ASlbar.*A15bar.*A33bar).*F3 + ... 

(Allbar.*A33bar - AlSbar.*A31bar)./(AlIbar.*A33bar.*ASSbar - 
*A3Sbar.*AS3bar - ASlbar.*A13bar.*ASSbar + ASlbar.*AlSbar.*AS3bar + 

*A13bar.*A3Sbar - ASIbar.*AlSbar.*A33bar).*FS; 
eta6_c = {A42bar.*A64bar - A44bar.*A62bar)./... 

{A22bar.*A44bar.*A66bar - A22bar.*A46bar.*A64bar - A42bar.*A24bar.*A66bar + 
*A26bar.*A64bar + A62bar.*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F2 - ... 

{A22bar.*A64bar - A24bar.*A62bar)./{A22bar.*A44bar.*A66bar - 
*A46bar.*A64bar - A42bar.*A24bar.*A66bar + A42bar.*A26bar.*A64bar + 

*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F4 + ... 

{A22bar.*A44bar - A24bar.*A42bar)./{A22bar.*A44bar.*A66bar - 
*A46bar.*A64bar - A42bar.*A24bar.*A66bar + A42bar.*A26bar.*A64bar + 

*A24bar.*A46bar - A62bar.*A26bar.*A44bar).*F6; 

% Absolute motions 


xil_point = etal_c 
xi2_point = eta2_c 
xi3_point = eta3_c 
xi3_prop = eta3_c 
xi3_mast = eta3_c 
xi3_bow = eta3_c 


+ etaS_c*z_point 
+ eta6_c*x_point 

- etaS_c*x_point 

- etaS_c*x_prop 

- etaS_c*x_mast 

- etaS_c*{L/2); 


- eta6_c*y_point; 

- eta4_c*z_point; 
+ eta4_c*y_point; 
+ eta4_c*y_prop; 

+ eta4_c*y_mast; 


forward at POI 
lateral at POI 
vertical at POI 
vertical at prop 
vertical at mast 
vertical at bow 


Relative motions/velocities/accelerations 


xi3_ 

vi3_ 

vi3_ 

xi3_ 

vi3_ 

ai3_ 

ai2_ 

vi3_ 

ai3_ 

xi3_ 

xi3_ 


bow_rel 

bow_rel 

.point 

.point_rel 

.point_rel 

.point 

.point 

.bow 

.bow 

.prop_rel 

.mast_rel 


= xi3_bow 


= xi3_point 


xi3_prop 

xi3_mast 


zeta; 

% 

relative 

vertical 

motion at 

bow 

xi3_bow_rel; 

% 

relative 

vertical 

velocity 

at bow 

xi3_point; 

% 

absolute 

vertical 

velocity 

at POI 

- zeta; 

% 

relative 

vertical 

motion at 

POI 

xi3_point_rel; 

% 

relative 

vertical 

velocity 

at POI 

\.'^ 2 ) .*xi3_point; 

% 

absolute 

vertical 

accel at 

POI 

\.'^ 2 ) .*xi2_point; 

% 

absolute 

lateral accel at POI 

xi3_bow; 

% 

absolute 

vertical 

velocity 

at bow 

\.^ 2 ) .*xi3_bow; 

% 

absolute 

vertical 

accel at 

bow 

• zeta; 

% 

relative 

vertical 

motion at 

prop 

• zeta; 

% 

relative 

vertical 

motion at 

mast 


% Begin random wave calculations 

% Long crested seas; Bretschneider spectral formulation 


% Sea spectrum 

S = {1.2S/4)*{ {omega_m'^4 . / omega .'^S)*{HS'^2)). *exp (-1.2S* {omega_m'^4 . / omega. ^4) ) ; 

Se = S./abs{{1-(2.0/g)*omega*V_calc*cos{beta_calc))); % Convert S(w) to S{we) 


% Define response spectra 


Seta3_c 

Seta4_c 

SetaS_c 

Sxil_point 

Sxi2_point 

Sxi3_point 

Sxi3_bow 

Sxi3_bow_rel 

Svi3_point 

Sai3_point 

Sai2_point 


( (abs (eta3_c) ) .'^ 2 ) . *Se; 

{(abs{eta4_c)) .^ 2 ) .*Se; 

{ (abs {etaS_c) ) . '^2) . *Se; 

{(abs{xil_point)).^2).*Se; 

{(abs{xi2_point)) .^ 2 ) .*Se; 

{(abs{xi3_point)) .^ 2 ) .*Se; 

{(abs{xi3_bow)) .^ 2 ) .*Se; 

{(abs{xi3_bow_rel)) .^ 2 ) .*Se; 
{(abs{vi3_point)) .^ 2 ) .*Se; 

{(abs{ai3_point)) .^ 2 ) .*Se; 

{ (abs {ai2_point) ) . "^2) . *Se; 
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Svi3_bow 

Svi3_bow_rel 

Sai3_bow 

Svi3_point_rel 

Sxi3_point_rel 

Sxi3_point_rel_2 

Sxi3_prop_rel 

S Xi3_prop_re1_2 

Sxi3_mast_rel 

Sxi3_mast_rel_2 


( (abs (vi3_bow) ) . ''2) . *Se; 

{{abs(vi3_bow_rel)) .^2) .*Se; 

((abs{ai3_bow)) .^2) . *Se; 

((abs{vi3_point_rel)).^2).*Se; 

{ (abs{xi3_point_rel)) .^2) .*Se; 

(omegae . "^2) . * { (abs {xi3_point_rel) ) . "^2) . *Se; 
{(abs{xi3_prop_rel)).^2).*Se; 

(omegae .^2) .*((abs(xi3_prop_rel)) .^2) .*Se; 

{{abs(xi3_mast_rel)) .^2) . *Se; 

(omegae .^2) .*((abs(xi3_mast_rel)) .^2) .*Se; 


% Initializations 


Seta3_c_i = 0; 
Seta4_c_i = 0; 
Seta5_c_i = 0; 
Sxil_point_i = 0; 
Sxi2_point_i = 0; 
Sxi3_point_i = 0; 
Sxi3_bow_i = 0; 
Sxi3_bow_rel_i = 0; 
Svi3_point_i = 0; 
Sai3_point_i = 0; 
Sai2_point_i = 0; 
Svi3_bow_i = 0; 
Svi3_bow_rel_i = 0; 
Sai3_bow_i = 0; 
Svi3_point_rel_i = 0; 
Sxi3_point_rel_i = 0; 
Sxi3_point_rel_2_i = 0; 
Sxi3_prop_rel_i = 0; 
Sxi3_prop_rel_2_i = 0; 
Sxi3_mast_rel_i = 0; 
Sxi3_mast_rel_2_i = 0; 


% Integral S (w) * | RAO | '^2 


for 1=2:1:filesize(2), 


Seta3_c_i = Seta3_c_i 

+ 

0.5* 

(Seta3_c(I) 

4 - 

Seta3_c(I-l)) 

* abs((omegae(I-l)-omegae(I))); 

Seta4_c_i = Seta4_c_i 

+ 

0.5* 

(Seta4_c(I) 

+ 

Seta4_c(I-l)) 

* abs((omegae(I-l)-omegae(I))); 

Seta5_c_i = Seta5_c_i 

+ 

0.5* 

(Seta5_c(I) 

+ 

Seta5_c(I-l)) 

* abs((omegae(I-l)-omegae(I))); 

Sxil_point_i = Sxil_point_i 

+ 

0.5* 

(Sxil_point(I) 

+ 

Sxil_point(I- 

1)) * abs((omegae(I-l)-omegae(I))) ; 

Sxi2_point_i = Sxi2_point_i 

+ 

0.5* 

(Sxi2_point (I) 

+ 

Sxi2_point(I- 

1)) * abs((omegae(I-l)-omegae(I))); 

Sxi3_point_i = Sxi3_point_i 

+ 

0.5* 

(Sxi3_point(I) 

+ 

Sxi3_point(I- 

1)) * abs((omegae(I-l)-omegae(I))); 

Sxi3_bow_i = Sxi3_bow_i 

+ 

0.5* 

(Sxi3_bow(I) 

+ 

Sxi3_bow(I-l)) 

* abs((omegae(I-l)-omegae(I))); 

Sxi3_bow_rel_i = Sxi3_bow_rel_i 

+ 

0.5* 

(Sxi3_bow_rel(I) 

+ 



Sxi3_bow_rel(I-l)) 


abs((omegae(I-l)-omegae(I))), 


Svi3_point_i = Svi3_point_i 

1)) * abs((omegae(I-l)-omegae(I))); 

Sai3_point_i = Sai3_point_i 

1)) * abs((omegae(I-l)-omegae(I))); 

Sai2_point_i = Sai2_point_i 

1)) * abs((omegae(I-l)-omegae(I))) ; 

Sai3_bow_i = Sai3_bow_i 

* abs((omegae(I-l)-omegae(I))); 

Svi3_bow_i = Svi3_bow_i 

* abs((omegae(I-l)-omegae(I))); 

Svi3_bow_rel_i = Svi3_bow_rel_i 
Svi3_bow_rel(I-l)) * abs((omegae(I-l)-omegae 

Svi3_point_rel_i = Svi3_point_rel_. 
Svi3_point_rel(I-l)) * abs((omegae(I-l)-omegae 
Sxi3_point_rel_i = Sxi3_point_rel_. 
Sxi3_point_rel(I-l)) * abs((omegae(I-l)-omegae 


+ 0.5*(Svi3_point(I) 


+ 0 , 

+ 0 , 

+ 0 , 

+ 0 , 

+ 0 , 
( 1 ) ) ) ; 
i + 0, 
( 1 ) ) ) ; 
i + 0, 
(I) ) ); 


5*(Sai3_point (I) 

5*(Sai2_point(I) 

5*{Sai3_bow{I) 

5*{Svi3_bow{I) 

5*{Svi3_bow_rel(I) 

5*{Svi3_point_rel(I) 
5*{Sxi3_point_rel(I) 


+ Svi3_point(I- 
+ Sai3_point(I- 
+ Sai2_point{I- 
+ Sai3_bow(I-l)) 
+ Svi3_bow{I-l)) 
+ 

+ 

+ 
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Sxi3_prop_rel_i = Sxi3_prop_rel_i + 0.5*(Sxi3_prop_rel(I) + 

Sxi3_prop_rel(I-l)) * abs((omegae(I-l)-omegae(I))); 

Sxi3_mast_rel_i = Sxi3_mast_rel_i + 0.5*(Sxi3_mast_rel(I) + 

Sxi3_mast_rel(I-l)) * abs((omegae(I-l)-omegae(I))); 

end 


% Second moment m_2 


for I=2:l:filesize(2), 

Sxi3_point_rel_2_i = Sxi3_point_rel_2_i + 0.5*(Sxi3_point_rel_2(I) + 
Sxi3_point_rel_2(I-l)) * abs((omegae(I-l)-omegae(I))); 

Sxi3_prop_rel_2_i = Sxi3_prop_rel_2_i + 0.5*(Sxi3_prop_rel_2(I) + 

Sxi3_prop_rel_2(I-l)) * abs((omegae(I-l)-omegae(I))); 

Sxi3_mast_rel_2_i = Sxi3_mast_rel_2_i + 0.5*(Sxi3_mast_rel_2(I) + 

Sxi3_mast_rel_2(I-l)) * abs((omegae(I-l)-omegae(I))); 

end 

% RMS values 


RMS_eta3_c 

RMS_eta4_c 

RMS_eta5_c 

RMS_xil_point 

RMS_xi2_point 

RMS_xi3_point 

RMS_xi3_bow 

RMS_vi3_point 

RMS_ai3_point 

RMS_ai2_point 

RMS_vi3_bow 

RMS_ai3_bow 

RM S_vi3_bow_re1 

RMS_xi3_bow_re1 

RMS_vi3_point_rel 

RMS_xi3_point_rel 

RMS_xi3_point_re1_2 

RMS_xi3_prop_re1 

RM S_xi3_prop_re1_2 

RMS_xi3_mast_rel 

RMS_xi3_mast_rel_2 


sqrt(abs(Seta3_c_i)); 

57.3*sqrt(abs(Seta4_c_i)); 

57.3*sqrt(abs(Seta5_c_i)); 
sqrt(abs(Sxil_point_i)); 
sqrt(abs(Sxi2_point_i)); 
sqrt(abs(Sxi3_point_i)); 
sqrt(abs(Sxi3_bow_i)); 
sqrt(abs(Svi3_point_i)); 
sqrt(abs(Sai3_point_i)); 
sqrt(abs(Sai2_point_i)); 
sqrt(abs(Svi3_bow_i)); 
sqrt(abs(Sai3_bow_i)); 
sqrt(abs(Svi3_bow_rel_i)); 
sqrt(abs(Sxi3_bow_rel_i)); 
sqrt(abs(Svi3_point_rel_i)); 
sqrt(abs(Sxi3_point_rel_i)); 
sqrt(abs(Sxi3_point_rel_2_i)); 
sqrt(abs(Sxi3_prop_rel_i)); 
sqrt(abs(Sxi3_prop_rel_2_i)); 
sqrt(abs(Sxi3_mast_rel_i)); 
sqrt(abs(Sxi3_mast_rel_2_i)); 


RMS_eta3_c_vector(ibeta,iV) 
RMS_eta4_c_vector(ibeta,iV) 
RMS_eta5_c_vector(ibeta,iV) 
RMS_xil_point_vector(ibeta,iV) 
RMS_xi2_point_vector(ibeta,iV) 
RMS_xi3_point_vector(ibeta,iV) 
RMS_xi3_bow_vector(ibeta,iV) 
RMS_vi3_point_vector(ibeta,iV) 
RMS_ai3_point_vector(ibeta,iV) 
RMS_ai2_point_vector(ibeta,iV) 
RMS_vi3_bow_vector(ibeta,iV) 
RMS_vi3_bow_rel_vector(ibeta,iV) 
RMS_xi3_bow_rel_vector(ibeta,iV) 
RMS_ai3_bow_vector(ibeta,iV) 
RMS_vi3_point_rel_vector(ibeta,iV) 
RMS_xi3_point_rel_vector(ibeta,iV) 
RMS_xi3_point_rel_2_vector(ibeta,iV) 
RMS_xi3_prop_rel_vector(ibeta,iV) 
RMS_xi3_prop_rel_2_vector(ibeta,iV) 
RMS_xi3_mast_rel_vector(ibeta,iV) 
RMS_xi3_mast_rel_2_vector(ibeta,iV) 


RMS_eta3_c; 
RMS_eta4_c; 
RMS_eta5_c; 
RMS_xil_point; 
RMS_xi2_point; 
RMS_xi3_point; 

RM S_xi3_bow; 

RMS_vi3_point; 

RMS_ai3_point; 

RMS_ai2_point; 

RMS_vi3_bow; 

RMS_vi3_bow_rel; 

RMS_xi3_bow_rel; 

RMS_ai3_bow; 

RMS_vi3_point_rel; 

RMS_xi3_point_rel; 

RMS_xi3_point_rel_2; 

RMS_xi3_prop_rel; 

RMS_xi3_prop_rel_2; 

RMS_xi3_mast_rel; 

RMS_xi3_mast_rel_2; 


N_P_prop(ibeta,iV)=3600*(1/sqrt(2*pi))*(RMS_xi3_prop_rel_2/RMS_xi3_prop_rel)*exp(- 
2*z_prop'^2/ (2*RMS_xi3_prop_rel) ) ; 

N_P_mast(ibeta,iV)=time_hrs*3600*(1/sqrt(2*pi))*(RMS_xi3_mast_rel_2/RMS_xi3_mast_rel)*exp 
(-2*z_mast'^2/ (2*RMS_xi3_mast_rel) ) ; 
end 


end 
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lL=larrLbdaoverL; 

temp_str = num2str(abs(str2num(beta_string)-180)); 

% Roll 
figure(1) 

[th,r]=meshgrid{{0:beta_incr:360)*pi/180,V_min:delta_V:V_max) ; 

[X,Y]=pol2cart{th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_eta4_c_vector,c_p_roll),colorbar 

title (['Roll (deg); Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' cm; SS ' 
S_S_string]) 

figure_string=[ ' graphs\roll\ ' ,Loading_Cond_string, , KG_string,,S_S_string, '. eraf ' ]; 

saveas(gcf,figure_string, ' eraf ' ); 

close 

% Pitch 

figure (2) 

[th, r]=meshgrid((0:beta_incr:360)*pi/180,V_min:delta_V:V_max); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_eta5_c_vector,c_p_pitch) , colorbar 

title ([' Pitch (deg); Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' cm; SS ' 
S_S_string]) 

figure_string=[ ' graphs\pitch\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '. eraf ' ]; 

saveas(gcf,figure_string, ' eraf ' ); 

close 

% Absolute Vertical Acceleration at POI 
figure(3) 

[th;. r ] =meshgrid ( (0 : beta_incr : 360) *p i/180; V_min : delta_V: V_max) ; 

[X;Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_ai3_point_vector/g,c_p_ver_acc) , colorbar 

title([ 'Absolute Vertical Acceleration (g) ; Loading = ' Loading_Cond_string ' %; KG = ' 
KG_string ' cm; SS ' S_S_string]) 

figure_string=[ ' graphs\ver_acc\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '.eraf'] 

saveas(gcf,figure_string, ' eraf ' ); 
close 

% Absolute Lateral Acceleration at POI 
figure(4) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,V_min:delta_V:V_max); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_ai2_point_vector/g,c_p_lat_acc) , colorbar 

title([ 'Absolute Lateral Acceleration (g) ; Loading = ' Loading_Cond_string ' %; KG = ' 
KG_string ' cm; SS ' S_S_string]) 

figure_string=[ ' graphs\lat_acc\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '.eraf'] 

saveas(gcf,figure_string, ' eraf ' ); 
close 

% Slam Velocity at bow 
figure(5) 

[th, r]=meshgrid((0:beta_incr:360)*pi/180,V_min:delta_V:V_max); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_vi3_bow_rel_vector,c_p_bow_vel),colorbar 
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title (['Slam Velocity (m/sec); Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' 
cm; SS ' S_S_string]) 

figure_string=[ ' graphs\slam_vel\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '.emf 

] ; 

saveas(gcf,figure_string, ' emf ' ); 
close 

% Slam Acceleration at bow 
figure(6) 

[th;. r ] =meshgrid { {0 : beta_incr : 360) *p i/180;. V_min : delta_V: V_max) ; 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_ai3_bow_vector/g,c_p_bow_acc),colorbar 

title (['Slam Acceleration (g) ; Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' 
cm; SS ' S_S_string]) 

figure_string=[ ' graphs\slam_acc\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '.emf 

] ; 

saveas(gcf,figure_string, ' emf ' ); 
close 

% Propeller Racing Events 
figure(7) 

[th, r]=meshgrid((0:beta_incr:360)*pi/180,V_min:delta_V:V_max); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',N_P_prop,c_p_prop),colorbar 

title (['Propeller Racing (/hr); Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' 
cm; SS ' S_S_string ' ; Prop. = ' z_prop_keel_string ' cm']) 

figure_string=[ ' graphs\race\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '_',z_prop 

_keel_string, '. emf ' ]; 

saveas(gcf,figure_string, ' emf ' ); 

close 

% Mast Submergence Events 
figure(8) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,V_min:delta_V:V_max); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',N_P_mast,c_p_mast_trip),colorbar 

title (['Mast Sub. (/trip); Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' cm; SS 
' S_S_string ' ; Range = ' range_string ' NM; Mast = ' z_mast_top_string ' cm' ]) 
figure_string=[ ' graphs\mast\ ' ,Loading_Cond_string,,KG_string,,S_S_string, '_',range_ 
string, '_' ,z_mast_top_string, '. emf ' ]; 
saveas(gcf,figure_string, ' emf ' ); 
close 

% Operating Envelopes (for area calcs and 01 processing) 

figure(9) 
colormap([0 0 0]) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,V_min:delta_V:V_max); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

contourf(X',Y',RMS_eta4_c_vector,[c_p_roll_max c_p_roll_max]) 
contourf(X',Y',RMS_eta5_c_vector,[c_p_pitch_max c_p_pitch_max]) 
contourf(X',Y',RMS_ai3_point_vector/g,[c_p_ver_acc_max c_p_ver_acc_max]) 
contourf(X',Y',RMS_ai2_point_vector/g,[c_p_lat_acc_max c_p_lat_acc_max]) 
contourf(X',Y',RMS_vi3_bow_rel_vector,[c_p_bow_vel_max c_p_bow_vel_max]) 
contourf(X',Y',RMS_ai3_bow_vector/g,[c_p_bow_acc_max c_p_bow_acc_max]) 
contourf(X',Y',N_P_prop,[c_p_prop_max c_p_prop_max]) 
contourf(X',Y',N_P_mast,[c_p_mast_trip_max c_p_mast_trip_max]) 

title (['01; Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' cm; SS ' S_S_string 
'; Range = ' range_string ' NM; Mast = ' z_mast_top_string ' cm; Prop. = ' 
z_prop_keel_string ' cm']) 
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figure_string= [' graphs\0I\ ', Loading_Cond_string,,KG_string,,S_S_string, '_',range_st 

ring,,z_mast_top_string,,z_prop_keel_string, '.bmp']; 

saveas(gcf,figure_string, 'bmp ' ); 

figl_data=imread(figure_string) ; 

suml_data=sum(sum(figl_data)); 

hold off 

y2_data=0; 

y3_data=l; 

x2_data=41904922; 

x3_data=59138077; 

x_data=suml_data; 

a_data={y3_data-y2_data)/{x3_data-x2_data); 

b_data={y2_data*x3_data-y3_data*x2_data)/{x3_data-x2_data); 
y_dat a=a_dat a * x_dat a+b_dat a; 
y_data=y_data+0.01; 
if y_data > 1 
y_data=l; 

end 

if y_data < 0 
y_data=0; 

end 

OI=100*y_data; 

OI_string = strcat{ ' 01= ' ,num2str(01,3), ' % ' ); 

%text(-17,10,0I_string) 

OI_vector(i_Sea_State)=100*y_data; 

clear figl_data 

close 

% Sea state loop ends 
end 

% Compute overall operability index 

001 = 100*sum((OI_vector/100).*(percent_probab/100)); 

OOI_string = strcat( ' 001= ' ,num2str(001,3) ) ; 
figure(10) 

colormap([0 0 1]),bar(0I_vector),xlabel( ' Sea State ') ,ylabel( ' Operability Index (%) ' ) 
axis([090 100]) 

title ([' Loading = ' Loading_Cond_string ' %; KG = ' KG_string ' cm; Range = ' 
range_string ' NM; Mast = ' z_mast_top_string ' cm; Prop. = ' z_prop_keel_string ' cm']) 
text(7,80,OOI_string) 

figure_string=[ ' graphs\00I\ ' ,Loading_Cond_string,,KG_string,,range_string, '_',z_mas 

t_top_string,,z_prop_keel_string, '.bmp']; 

saveas(gcf,figure_string, 'bmp ' ); 

close 

clear all 
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APPENDIX B. SHIPMO SAMPLE INPUT EILE 


ASCC 
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APPENDIX C. NONCRITICAL OUTPUT RESULTS 



Figure 60. Slam acceleration at 80% midbox load. 



Figure 61 . Slam acceleration at 50% midbox load. 
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Figure 62. Lateral acceleration at 80% midbox load. 
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Figure 64. 


Figure 65. 
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Figure 66. 


Figure 67. 
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